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Abstract. We obtain threshold conditions for eradication and permanence 
of the disease for a nonautonomous SIRQ model with time-dependent inci- 
dence rates, given by functions of all compartments in some general family. 
The threshold conditions are given by some numbers that play the role of the 
basic reproduction number. Additionally, we obtain simple threshold condi- 
tions in the asymptotically autonomous and periodic settings and show that 
our thresholds coincide with the basic reproduction number in some common 
autonomous examples. 



1. Introduction 

Quarantine is one of the possible measures for reducing the transmission of dis- 
eases and has been used over the centuries to reduce the transmission of human 
diseases such as leprosy, plague, typhus, cholera, yellow fever and tuberculosis. 

The concept of quarantine dates back to the fourteenth century and is related 
to plague. In fact, in 1377, the Rector of the seaport of the old town of Ragusa 
(modern Dubrovnik) officially issued a thirty-day isolation period: ships coming 
from infected or suspected to be infected sites were to stay at anchor for thirty 
days before docking. A similar measure was adopted for land travelers but the 
period of time was enlarged to forty days. It is from the Italian word for forty - 
quaranta - that comes the term quarantine [9] . 

Isolation is still nowadays an effective way of controlling infectious diseases. For 
some milder infectious diseases, quarantined individuals could be individuals that 
choose to stay home from school or work while for more severe infectious diseases 
people may be forced into quarantine. 

To study the effect of quarantine in the spread of a disease, we consider a SIRQ 
model. Thus, we divide the population into four compartments: the S compartment 
corresponding to not infected individuals that are susceptible to the disease, the / 
compartment corresponding to individuals that are infected and not isolated, the Q 
compartment corresponding to isolated individuals (Q stands for quarantine) and 
the R compartment including the recovered and immune individuals. We assume 
in this paper that the quarantined individuals are perfectly isolated from the others 
so that they do not infect the susceptibles and that the infection confers permanent 
immunity upon recovery. 

A major concern when studying epidemic models is to understand the asymptotic 
behavior of the compartments considered, particularly the infective compartment. 
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Frequently, in autonomous epidemic models, one of two situations occur: there is a 
disease-free equilibrium that is the unique equilibrium of the system or the disease- 
free equilibrium coexists with an endemic equilibrium. The basic reproduction 
number is a fundamental tool to obtain the behavior for a given set of parameters: 
if the basic reproductive number is less than one the disease dies out and if this 
number exceeds one the disease remains in the population. 

Several authors studied autonomous epidemic models with a quarantine class. 
Namely, Feng and Thieme [5] proposed a SIQR model for childhood diseases and 
showed that isolation can be responsible for the existence of self-sustained oscilla- 
tions. Their model includes a modified incidence function - the quarantine-adjusted 
incidence - where the number of per capita contacts between susceptibles and infec- 
tives is divided by the number of non isolated individuals. Some dynamical aspects 
of this SIQR model were studied by Wu and Feng [12]. Namely, these authors 
studied this model through unfolding analysis of a normal form derived from the 
model and found Hopf and homoclinic bifurcations associated to this unfolding. A 
generalized version of Feng and Thieme's model, allowing infected individuals to 
recover without passing through quarantine and also disease related deaths, was 
proposed by Hethcote, Zhien and Shengbing in [5], where this model was compared 
to SIQR models with different incidence rates. A similar analysis is also under- 
taken in that paper for SIQS models with different incidence rates. The authors 
concluded in that work that the SIQR model with quarantine-adjusted incidence 
seems consistent with the sustained oscillations that are observed in real disease 
incidences but do not match every aspect of the observed data. 

Other types of epidemic models with quarantine were studied in the literature. 
A SEIQR model was considered by Gerberry and Milner in [3] where the authors 
studied the dynamics and used historical data to discuss the adequacy of the model 
in the context of childhood diseases. The global dynamic behavior of a different 
SEIQS epidemic model, with a nonlinear incidence rate of the form j3SI, was ana- 
lyzed by Yi, Zhao and Zhang in |17j . In another direction, Arino, Jordan and van 
den Driessche [I] studied the effect of quarantine in a general model of a disease 
that can be transmitted between different species and multiple patches (though 
in their setting quarantine is present only in the form of travel restriction between 
patches). A model with twelve different classes was proposed by Safi and Gumel [16] 
to study the impact of quarantine of latent cases, isolation of symptomatic cases 
and of an imperfect vaccine in the control of the spread of an infectious disease in 
the population. The same authors studied a delayed model including a quarantined 
class [15] . 

We emphasize that all the above models are autonomous models. However, it is 
well-known that fluctuations are very common in disease transmission. For instance, 
the opening and closing of schools tend to make the contact rates vary seasonally in 
childhood diseases [5] and periodic changes in birth rates of the population are also 
frequent. It is therefore natural to consider time-dependent parameters in epidemic 
models and define threshold conditions for extinction and persistence of the disease. 

Our objective in this paper is to consider a family nonautonomous SIQR models 
with a large family of contact rates given by general functions that are, in general, 
time-dependent and far from bilinear, and obtain numbers that play the role of 
the basic reproductive number in this setting. Related results for nonautonomous 
SIRVS models were obtained in [IT] [13] , thought the incidence rates considered in 
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those papers are linear in the infectives while here we assume that the incidence 
rates are, in general, nonlinear functions of all compartments. 

The structure of the paper is as follows: in section[5]we introduce the model and 
the main definitions; in section [3] we state our main Theorem; in the next sections 
we look at some particular cases of our model, namely, in section|4l we particularize 
for the autonomous model, in section [SJ we see that the asymptotically autonomous 
model has the same thresholds as the limiting autonomous model and in section [6] 
we obtain thresholds for the periodic model with constant death and recruitment 
rates; finally, in section [7] we prove some auxiliary results and the main theorem. 

2. Generalized nonautonomous SIQR model 

We consider the following nonautonomous and general incidence rate SIRQ 
model that generalizes the three autonomous SIRQ models in [5], 

' S' = A(i) - (p(t, S, R, Q, I) - d(t)S 

I' = p(t, S,R,Q,I)- [ 7 (t) + a{t)+ d{t) + a x (*)] I . . 

a(t)I - [d{t) + a 2 (t) + e(t)]Q 
R' = j{t)I + e(t)Q - d{t)R 

where S, R, Q and I correspond respectively to the susceptible, recovered, quaran- 
tine and infective compartments; A(t) is the recruitment rate of susceptible (births 
and immigration); tp(t, S, R, I) is the incidence (into the infective class); d(t) is the 
per capita natural mortality rate; 7(f) is the recovery rate; a(t) is the removal rate 
from the infectives; otx(t) is the disease related death in the infectives; ot2(t) is the 
disease related death in the quarantine class; and s(t) is the removal rate from the 
quarantine class. Furthermore we assume that A, d, 7, a, ax, eti and e are continuous 
bounded and nonnegative functions on Rj and that 

inf d(t) > and A = supA(i)>0. 



t>0 

For each 6 £ Rj define the set 

A s = {(xx,X2,x 3 ,X4,) el 4 : <x t < A/d + S}. 

Assume further that the function ip : (Rj) 5 — > R is continuous and nonnegative 
and that there are K > and M > such that, given < 6 < M, we have 

<p(t,0,y,w,z)=0, (2) 

\ip(t,xx,y,w,z) - ip(t,x 2 ,y,w,z)\ < K\x x - x 2 \z, (3) 

ip(t,x,0,0,r) (p(t,x,y,w,z) . ip(t, x, 0,0, t) 
sup > > mi , (4) 

0<r<\/d+S T Z 0<r<X/d+S T 

for all (t, 0, y, w, z), (t,x,y,w,z), (t, xx, y, w, z), (t, xx, y, w, z) e R(j~ x A$. Note 
that, by ([2]) and ([3]), the above infimums and supremums are finite. Note also that, 
for contact rates that are not dependent on y and w, ([4]) is always verified. 
Now, we need to make some definitions. We say that 

i) the infectives / are permanent in the solution (S(t),I(i), R(t),V(t)) of ([1]) if 
there exist constants < mi < m 2 such that 

mi < liminfl(t) < limsupl(t) < m,2', 
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ii) the infectives I go to extinction in the solution (S (t) , I (t) , R(t) , V(t)) of §T§ if 
lim I(t) = 0; 

t— >oo 

iii) the infectives / are permanent in system ([1]) if there exist constants < mi < 
?7i2 such that 

mi < liminf I(t) < lim sup /(i) < m2, 

t-»-oo t->oo 

for every solution (S(t),I{t),R(t), V(t)) of ([I]) with Sb(t), h{t), Ro(t),V (t) > 
(note that we require that mi and m-2 are independent of the given solution 
with positive initial condition in /); 

iv) the infectives / go to extinction in system ([1} if lim I(t) = for all solution 

t— >oo 

Of ©. 

3. Statment of the results 
Consider the auxiliary equation 

x > = A(t) - d(t)x, (5) 
and, for each solution x(t) of ([5]) with x(0) > 0, the functions 

b I , s (t,x(t))= inf ^ °' °' T) - ( 7 (t) + <x(t) + d(t) + ai (t)) 

0<t<^+<5 T 

and 

&s,*(t,a;(i)) = sup ^ X ®' °' °' r) - ( 7 (t) + tr(t) + d(t) + ai(t)) ■ 

0<T<|+5 T 

Finally, define the numbers 

pt+X 

r p (A) = lim liminf / bi g(s,x(s)) ds, 

<5— s-0 t— >CX3 J t 

/■t+A 

r e i (A) = lim lim sup / b$ s{s, x(s)) ds, 

<5->0 t^oo J t 

1 I* 

r e ,2 = lim limsup — / bs,s{s, x(s)) ds. 

5->0 t->oo t Jq 

and 

JJp(A) = e r ^ x \ R e ,x(X) = e r "- lW and R e<2 = e r ^ 2 . 
Contrarily to what one could expect, the next technical lemma that will be proved 
in Section [7] shows that the numbers above do not depend on the particular solution 
x(t) of §5§ with x(0) > 0. The numbers R P (X), R e ,i(X) andi? e ^ play the role of 
the basic reproduction number (defined for autonomous systems and also called 
quarantine number) in the non-autonomous setting. 
We have the following result. 

Lemma 1. We have the following: 

1. Let e,ex,e 2 > and < 5 < M. If 

a, b, a ± ex, b ± £2 S ]0, A/d+ S[ and a — b<£, 

then 

\b It5 (t, a ± ex) - bi, s (t, b ± e 2 )\ < (e + £1 + e 2 )K (6) 

and 

\b s ,s(t,a±ei) -b s ,s(t,b±e 2 )\ < (e + e 1 +£ 2 )K. (7) 
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2. The numbers R p {\), R e s{\) and i? e ,2 are independent of the particular so- 
lution x(t) with x(0) > of ([5]). 

We now state our main theorem on the permanence and extinction of the infec- 
tives in system ([T]). 

Theorem 1. We have the following for the system ([TJ. 

1. If there is a constant A > such that R P (X) > 1 then the infectives I are 
permanent. 

2. If there is a constant X such that R e ^i(X) < 1 or if R e .2 < 1 then the infectives 
I go to extinction and any disease free solution (S(t), 0, R(t), Q(t)) is globally 
attractive. 

4. Autonomous model 

In this section we are going to consider the autonomous setting. That is, we are 
going to assume in system ([lj that A, d, 7, 8, ax, e and a 2 are constant functions, 
that A, d > 0, and that ip is independent of t. We obtain the autonomous system 

f S' = A- (fo(S, R, Q, I) - dS 
V = <p (S, R,Q,I)-[j + 5 + d + ai }I 
Q' = 5I-[d + a 2 +e]Q 
^R' = 7/ + eQ - dR 

where ipo '■ (Ito~) 4 — ^ IR is continuous, nonnegative and satisfies 

(£o(0, y,w,z) = 0, \(p (xi,y,w,z) - ip (x 2 ,y,w, z)\ < K\x\ -x 2 \z, 

and 

(Po{x,0,0,t) 



(8) 



sup 

a<r<\/d+s 



(Po(x,0,0,t) ip (x,y,w,z) 



inf 

0<T<\/d+S 



for all (0, y, w, z), (xx, y, w, z), [x 2 , y, w, z) G As. In this setting we have that the 
auxiliary equation ([5]) admits the constant solution x(t) = A/d. Thus we obtain 



r P (A) = 



. y> o (A/d,0,0,T) 
mi (7 + + d + ax) 

0<r<4 T 



^e,l(A) 



and 



r e ,2(A) 



0<T<i 



sup 

0<r<^ 



A, R P (X) = e f >< A >, 
A, i? 6i i(A) 

MmOAr) _ h + s + d + aih Re2(x) = ercMX} 



y>o(A/rf,0,0,r) 
sup (7 + + d + ax) 



= e r.,!(A) 



It is now easy to establish a result that is a version of Theorem Q] in the particular 
case of autonomous systems. Define 



R 



AUT 

■p 



= inf 



¥>o(A/d,0,0,T) 



o<r<£ (7 + 5 + d + ax)n 



and R£ UT = sup 



0<T< 



y>o(A/rf,0,0,T) 
a (7 + S + d + ai)r 



For any A > 0, we have R P (X) > 1 if and only if R p UT > 1 and i? ej i(A) < 1 if and 
only if Rf UT < 1 and this implies the following result. 

Theorem 2. We have the following for the autonomous system ([5]). 
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1. If Rp UT > 1 then the infectives I are permanent; 

2. If R^ T < 1 then the infectives I go to extinction and any disease free 
solution is globally attractive. 

Note that, if the incidence rate has the particular form 

<Po(S, R, Q, I) = ip(S, R, Q)g(I)I, (9) 

where g is a continuous, bounded and nonnegative function, then 

ip(A/d,0,0)A inf g(r) ip(A/d,Q,0)A sup g{r) 

Rt UT = 77 ; °. <T< * and R^ UT — ° <T<4 



d(j + S + d + ai ) " e d(7 + (5 + <i + Q;i) 

and we conclude that the infectives are permanent (resp. go to extinction) if 

mf ff(r) > I,,.,, n n x (resp. sup g(r) < 



0<r<f^ ' A^(A/d,0,0) v r o<r<A" v y A^(A/rf, 0, 0) 7 
Naturally, if 

d(7 + S + d + a x ) /r« » ; n\ 

Vmo,o) £g(|0 ' A/dD 

we get no information about the asymptotic behavior of the infectives. 

Assuming now that g is constant, we have Rp UT — R^ UT and we can re- 
cover the results obtained in [5]. In fact, the autonomous SIQR model with 
simple incidence (ipo(S, R,Q, I) = /3SI), the autonomous SIQR model with stan- 
dard incidence ((po(S,R,Q,I) = f3SI/(S + I + Q)) and the autonomous SIQR 
model with quarantine-adjusted incidence {fo{S, R, Q, I) — 0SI / (S + I + Q + R)) 
are all in the conditions we have assumed and our threshold values for those 
model are Rp UT = Rf UT = /3A/[d(7 + S + d + a±)] for simple incidence and 
Rp UT = Rf UT = + 5 + d + a{) for standard and quarantine-adjusted inci- 
dence. For these models, the numbers Rp UT = R^ UT coincide with the quarantine 
reproduction number obtained in [5]. Note that, to put the referred autonomous 
models considered in [5] in our context, we need to have the standard incidence and 
the quarantine-adjusted incidence continuous in A$. This constitutes no problem 
since the functions (p associated to these incidence rates can be extended to a 
continuous function in A^, setting ifo(0, i?, 0, 0) = for the standard incidence and 
c/?o(0, 0, 0, 0) = for the quarantine-adjusted incidence. 

In some situations, we may be able to obtain information about extinction but 
not about permanence or vice versa. To illustrate this fact we consider the family 
of contact rates obtained by setting ip(S, Q, R) = S and g(I) = 1 ^ Iq with p, q > 0, 
a > 0, in (|9|). This family of contact rates was considered for instance in [T4l [6], 
Our threshold conditions above give no relevant information about the permanence 
of the disease but it is easy to see that if 

(A/d) p - 2 - [1 + a(A/d) 9 " 4 ] (7 + S + d + ax) < 

then the disease goes to extinction. 

5. Asymptotically autonomous model 



In this section we are going to consider the asymptotically autonomous SIRQ 
model. That is, additionally to the assumptions on Theorem Q] we are going to 
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assume for system {!]) that there is a continuous function tpo such that 
lim <p(t, S, R, Q, I) = <p {S, R, Q, I) 

t— »oc 

for each (S, R, Q, I) € R 4 and that the time-dependent parameters are asymptot- 
ically constant: /3(t) -» (3, A(t) -» A, d(t) ->■ d, j(t) -> 7, a(t) -> <5, ai(i) -> ai, 
«2(i) — ► «2 and e(t) — »■ £ as < — > +00. Denoting by F(t, S 1 , i?, Q, J) the right hand 
side of ([1]) and by Fo(5, R, Q, I) the right hand side of the limiting system, that is 
of JSJ, we also need to assume that 

lim F(t,S,R,Q,I)=F (S,R,Q,I), 

t— »+oo 

with uniform convergence on every compact set of (H^) 4 and we will also assume 
that (S,R,Q,I) h-> F(t,S,R,Q,I) and (S,R,Q,I) h-> F (S,R,Q,I) are locally 
Lipschitz functions. 

There is a general setting that will allow us to study this case. 

Let / : R x R n —> R and fo : R n R be continuous and locally Lipschitz in 
R™. Assume also that the nonautonomous system 

x' = f(t,x) (10) 

is asymptotically autonomous with limit equation 

x 1 = fo(x), (11) 

that is, assume that f(t,x) — > g(x) as t — » +00 with uniform convergence in every 
compact set of R™. The following theorem is a particular case of a result established 
in [TU] (for related results and applications see for example 0[7])- 

Theorem 3. Let $(i, to, xq) and ip(t,to,yo) be solutions of (jTUJ) and (jTTJ) respec- 
tively. Suppose that e <E R™ is a locally stable equilibrium point of (jTTJ) with attrac- 
tive region 

W(e) = !y€R n : lim <p(t,to,v) = e\ 

and that W$ PI W(e) ^ 0, where W$ denotes the omega limit of <&(i, to, xq). Then 
lim $(i, to, a?o) = e. 

Since (R + ) 4 is the attractive region for any solution of system ((5J with initial 
condition in (R + ) 4 and the omega limit of every orbit of the asymptotically au- 
tonomous system with I (to) > is contained in (R + ) 4 , we can use Theorem [3] to 
obtain the following result. 

Theorem 4. We have the following for the asymptotically autonomous systems 
above. 

1. If Rp UT > 1 then the infectives I are permanent; 

2. If Rf UT < 1 then the infectives I go to extinction and any disease free 
solution is globally attractive. 

6. Periodic model with constant natural death and recruitment 

In this section we are going to consider a periodic SIRQ model. Additionally to 
our assumptions on the function tp and the parameter functions in Theorem [U we 
are going to assume in system {1} that X(t) — A, d(t) — d are constant functions, 
that there is a T > such that <p(t, S, R, Q, I) = <p(t + T, S, R, Q, I) and that the 



8 



CESAR M. SILVA 



remaining time-dependent parameter functions are periodic functions with period 
T. We have in this case the constant solution x(t) = A/d and therefore 

r p (A) - Jim liminf f +T sup <P&AMlMll} _ ( 7 ( s ) + s(s) + d + a x (s)) ds 

0<r<±+5 T 



5->0 t-v+oo 



<p(s, A/d,0,0,r) , ,_ ? , _ NrT1 
sup ^ v / ' ' ' ' ds-ij + S + d + a^T 

o 0<t<A T 

and similarly 

r el (X)=[ T inf < fa-( ? + g + d + ai )r, 

Jo o<^<f T 

where / denotes the average of / in the interval [0,T]: f — ^ Jq f(s) ds. Define 
and 



o <r<A r(7 + 5 + d + ai) 



Theorem 5. We have the following for the periodic system with constant recruit- 
ment and death rates. 

1. If R PER > 1 then the infectives I are permanent; 

2. If R PER < 1 then the infectives I go to extinction and any disease free 
solution is globally attractive. 

If we assume in the periodic model above that the incidence rate has the partic- 
ular form ip{t, S, R, Q, I) = (3{t)SI with fi(t + T) = f3(t), then we obtain 

R PER {\) = K ER W = Tr T^l T- ( 12 ) 

di-f + S + d + m) 

and if we assume in that model that <p(t, S, R, Q, I) = f3(t)SI/(S + R + Q + I) or 
that (p(t, S, R, Q, I) = [3(t)SI/(S + R + I) with (3{t + T) = 0(t) we get 

R PER (X) = R PER (X) = —j^j— , (13) 
7 + o + a + ai 

and, assuming that the parameter functions are all constant, by (|12p and (|13p we 
obtain again the thresholds for the autonomous system with simple, standard and 
quarantine-adjusted incidence. 

7. Proof of Theorem [T] 

We will divide the proof in several steps. Namely, first we establish some simple 
results about the solutions of (JTJ, then we analyze the non- homogeneous linear 
system ([5]) and finally we prove the theorem. 

Proposition 1. We have the following: 
i) all solutions (S(t), I(t), Q{t), R(t)) of (JTJ) with nonnegative initial conditions 

are nonnegative for all t > 0; 
ii) all solutions (S(t), I(t),Q(t), R(t)) of ([lj with positive initial conditions are 
positive for all t > 0; 
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Hi) If (S (t) , I (t) , Q (t) , R(t)) is a solution of ([T]) with nonnegative initial conditions 
then 

timsupS(t) + I(t) + Q(t) + R(i) < A/d. 

t— > + 00 

Proof. Properties El and|H| are easy to prove. In fact, since t M- A(t) and t h-> /i(t) 
are bounded, setting U(t) = S(t) + I(t) + R(t) + V(t) and adding the equations 
in (JT|) we obtain for nonnegative initial conditions, 

U' = A{t) - d(t)U - a x (t)I - a 2 (t)Q < A — dll, 
where A = sup te]R + A(t) and d — inf tgE + d(t). Moreover, since the general solu- 
tion of equation V' = A — dV is V(t) = C e~ dt +A/d, we obtain by comparison 
property [m]). □ 

Proposition [T] shows that, for every e > 0, every solution (S (t) , I (t) , Q (t) , R(t)) 
of our system stays in the region 

B s = {(x,y,w,z) e (K+) 4 : x + y + w + z < A/d + S} C A s 

for every t € sufficiently large. 

For to > 0, the solution of ([5]) with initial condition x(to) — xq > is given by 

x(t) = e~ d(s) ds x + /* e" /« d ^ ds A(u) du 

J to 

and this allows us to establish easily the following proposition. 

Proposition 2. FFe /ia«e t/ie following: 
i) Given to > 0, all solutions x(t) of equation ([5]) with initial condition a; (to) > 

are nonnegative for all t > 0; 
n) Given to > 0, aZ/ solutions x(t) of equation ([5]) wit/i initial condition x(to) > 

are positive for all t > 0; 
to) Each fixed solution x(t) of ([5]) with initial condition x(to) > is bounded and 

globally uniformly attractive on [0, +oo[; 
iv) if x(t) is a solution of ([5]) and x(t) is a solution of the system 

x' = A(t) -d(t)x + f(t) 

with x(to) = x(to) then 

sup \x(t) - x(t)\ < 1/d sup |/(t)|. 

t>t t>t 

7.1. Proof of Lemma [H Assume that £,£i,£2 > 0, < 6 < M, a, b, a ± £i, 

6 ± e 2 €]0, A/d + <5[ and that a - 6 < s. We have, by 10, 

|y>(t,a±£i,0,0,T) - y(t, 6 ±e 2 ,0,0,r)| < K\a - b ± e x ± e 2 \r < Kr(e + e\ + £ 2 ). 
Therefore, 

¥>(t,a±£i,0,0,T) <y9(t,6±£ 2 ,0,0,T) 

A (£ + £j + £ 2 ) < 

T T (14) 

„ y>(t, ffl ± ei ,0,0,T) V ^ 

< h AT(£ + £i +£ 2 ), 

T 

and thus, taking the infimum and the supremum in r g]0, A/d + S[ in equation (| 14[) 
and adding and subtracting 7(t) + <r(t) + d(t) + ai(t), we get © and ([7]). 

We will now show that in fact r p (A), r e ,i(A) and r ej2 (A) are independent of the 
particular solution x(t) of ([5]) with x(0) > 0. In fact, by[m]) in Proposition [21 for 
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every e > and every solution x\(t) of (O with xi(0) > there is a T £ > such 
that \x(t) —Xx(t)\ < e for every t > T e . Let <5 G]0, M[ and choose es > such that, 
for e < £5 and t > T e , we have x(t), x\(t), x(t) ±e 6 ]0, A/d + 6[. Leting a — x(t), 
b = xi(t) and ei,£2 = we obtain by ((6|) 

/f + A rt + \ rt + \ 

bi^s(s,xi(s)) ds — XKe < / bi t $(s, x(s)) ds < / bi t $(s, x\ (s)) ds + XKe, 

for every t > T e . We conclude that, for every < e < es, 

rt+\ /-t+A 

liminf / hi s (s, x(s j) ds — liminf / 6/ <s(s, xi(s)) e?s < XKe, 

t^ + oo J t t^+oo J t 

and thus 

/-t+A /.t+A 

liminf / bj g (s, x(s)) ds = liminf / bi s(s, xi(s)) ds. 
t— >+oo y. t— ^+oo 

Letting <5 — ^ we conclude that r p (X) (and thus i? p (A)) is independent of the 
chosen solution. Using ((7]), the same reasoning shows that r e ,i(A) and r ej2 (and 
thus i? ej i(A) and R e ,2) are also independent of the particular solution. This proves 
the lemma. 

7.2. Proof of Theorem [U Assume that R P (X) > 1 (and thus r p (X) > 0) and let 
(S(t),I(t),R{t),Q(t)) be any solution of © with S(T ) > 0, I(T ) > 0, R(T ) > 
and Q(T ) > for some T > 0. 

Let x*(t) be some fixed solution of ([5]) with x*(0) > 0. Then, by[u| in Proposi- 
tion [2] we have x*(Tq) > 0. Since r p (X) > 0, using (|6|) we conclude that there are 
constants (5i e]0, M[, AT > and Ti > such that 

rt+A 

&/ :(5l (s,x*(s)-6>)ds> A (15) 



for all t > T\ and 6> > sufficiently small, say for < 9 < e. By [m]) and [rv]) in 
Proposition [T] we may assume that for t > To we have {S(t), I(t), R(t), Q(t)) £ . 
Define 

£o = min {^Ty s } and £l = 2^(A% d + ( 5 1 )- (16) 

We will show that 

limsupi(t) > e x . (17) 

t— >+oc 

Assume by contradiction that (|T7f is not true. Then there exists T 2 > satisfying 
I(t) < e\ for all t > T 2 . Consider the auxiliary equation 

x' = A(t) - d(t)x - K(A/d + 5i)ei. (18) 

Let x(t) be the solution of ([3]) with x(Tq) = S(To) and let x(t) be the solution 
of HHD with x{T Q ) = S(T ). By Ev]) in Proposition [1 and (HU), for all t > T , we 
have 

\x(t)-x(t)\<^/d + Si)ei = ^. (19) 

According to[mJ in Proposition [21 x*(t) is globally uniformly attractive on Rq . 
Thus there exists T3 > such that, for all t > T3, we have 

\x(t)-x*(t)\<^. (20) 
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From our assumptions about the function ip, we have 

tp(t, S(t), R(t),Q(t), I(t)) < KS{t)I{t) < K(A/d + Si)ei, 

for all t > max{To,T2}, and thus, by the first equation in ([T]), 

S' > k{t)-d{t)S -K{k/d + 5i)ei. (21) 

for all t > max{T ,T 2 }. Comparing (JTSJl and ([21} we conclude that S(t) > x(t) for 
all t > max{r ,T 2 }. Write T 4 = m&x{T 0l T ll T 2 ,T 3 }. By (JTSJ and ((201) we have, 
for all i > T 4 , 

S(t) > x(t) > x(t) - ^ > x*(t) - e„. (22) 
By and © in Lemma Q] we get, for all t > T 3 , 

b I (t,S(t))>b I (t,x*(t)-e )-2Ke . (23) 
According to (TTB1) we have £o < £ an d by (fT5)) we get 

^ 6/ A (r,x*(r) - £o )dr > Q(t-T 3 )-l^iV. (24) 

By the second equation in ([I) we obtain 
|>(t,S(t),fl(t),Q(t),I(t)) 



(7(*)+^(*) + d(*)+ai(*)) 



J(t) 



Z(i) 

and thus, integrating from T3 to t and using ((4]), ([23]) . (f24|) and (|16l) . we have 
I{t) = I(T 3 ) e^3 •P( s ' s ( s )' R ^),Q(B),i(s))/i(. s )-( J ( s )+6(s)+d( s )+ ai (s)) ds 

> 7(T 3 ) b '- Sl ( T ' s ( s ^ ds 

> I(Ts) B^3 b '.«l( T > x *( s ) _e °) ds-2Ks (t-T 3 ) 

> I(T 3 )c ( ^- 2K£a)(t - T3) - N 
>I(T 3 )e^ (t ~ T ^- N , 

and we conclude that I(t) — > +00. This contradicts the assumption that I(t) < E\ 
for all t >T 2 . From this we conclude that (fT7| holds. 

Next we will prove that for some constant £ > we have in fact 

liminfl(t) > £ (25) 

t— >-\-oo 

for every solution (S(t),I(t),R(t),Q(t)) with S(T Q ) > 0, 7(T ) > 0, i?(T ) > 
and Q(T ) > 0. By (JTSJ, there is a A x > A such that for all £ > Ai, t > T x and 
< 9 < e, we have 

b IiSl {s,x*(s)-9)ds> N. (26) 



We proceed by contradiction. Assume that (l2~5j) does not hold. Then there 
exists a sequence of initial values (x n ) n ^, with aj„ = (S n , R n , Q n , I n ) with S n > 0, 
i?„ > 0, Q n > and /„ > such that 

liminf lit, x n ) < 

t->+oo n 
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where I(t,x n ) denotes the solution of (JTJ with initial conditions S(T\) = S n , 
J(Ti) = I n , R(Ti) = R n and Q(Ti) = Q n . By (JTZJ) , given n 6 IN, there are 
two sequences (t ni k)kew and (s n ,fe)feeK with 

T\ < S n ,l < t n ,l < Sn,2 < *n,2 < • • • < S n ,fc < in.fc < ' ' ' 

and lim s n /. = +00, such that 

I(s fc ,„,X„) = — , I(tk,n,Xn) = -\ (27) 

and 

^ < I(t,x n ) < — , for all i 6]s ni fc,*n,fc[- (28) 
n n 

By the second equation in ([TJ we have 

J'(t, x n ) = I(t, x n )[(<p(t, S(t, x n ), R(t, x n ), Q(t, x n ),I(t, x n ))/I(t, x n )- 
- ( 7 (t) + S(t) + d(t) + ai(t))} 

> -( 7 (i) + cr(t) + d(t) + ai(t))/(t, i„) 

> - 7o /(t, £„), 

where 70 = sup s>0 { 7 (s) + <5(s) + d(s) + a±(s)}. Therefore we obtain 

r fc,™ j (r, x„) , 

— c(T > - 7 o(tfe,n - Sfe,n) 

»*,„ A 7 "^™) 

and thus I{t k>n ,x n ) > I{s k , n ,x n ) e -»(**.»-»*.»). By (27} we get ± > e -^( tfc --^.") 
and therefore we have 

logn 

tk,n-Sk.n> >• +OO (29) 

7o 

asn-> +00. 

According to ([2"5jl. for all i s]sfe, n , ijfe,n[) we have I(t,x n ) < ei/n < s\ and we 
conclude that ip(t, S, R,Q,I) < K(A/d + 5 1 )e 1 for all i e]sfe, n , i fcj „[. Thus by (JTJ) 
we obtain 

S" > A(t) - d(t)5 - if (A/d + <5i)ei 

for all t 6]s fc , n ,tfc, n [. 

By comparison we have S(t,x n ) > x{t) for all t &]sk, n ,tk, n [, where x is the 
solution of (US]) with x(sk, n ) — S(sk,n, x n)- Using (fH?|) we obtain, for all t 6 

|z(i)-z(i)|<|, (30) 

where x(i) is a solution of ([5]) with x(sk, n ) — S'(sfe ! „, x„). Since x*(t) is globally 
uniformly attractive, there exists T* > 0, independent of n and /c, such that 

\x(t)-x*(t)\< F f, (31) 

for all t > Sk, n + T* . By (|29|) we can choose B > such that — s n ,fc > Ai + T* 
for all n > B. Given n > B, by (|2^|) . (j3"0")) and (|31l) and by the second equation 
in (P) we get 

^ £ 1 /„*.*' n x-r« b I,Si (T,x*(T)-e )dT _ £1 
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This leads to a contradiction and establishes that liminf I(t) > £ > 0. Thus, [I] in 

Theorem Q] is established. 

Assume now that i? e ,i(A) < for some A > or that R e ^ < 0. Then we can 
choose Si G [0, M], N > and T > such that 

t+A 

& s , 4l (s,a:*(s) + 0)ds < -JV (32) 



or 

rt 



b s , Sl {s,x*(s)+9)ds<-N (33) 



1 

7 

for all t > Tq and > sufficiently small. By [m]) in Proposition [T] we can assume 
that (S(t),I{t),R(t), Q{t)) e A 5l for i > T . Given a solution of © with S(T ) > 0, 
fl(T ) > 0, Q(T ) > and I(T ) > we have 

S'(t) <A(t)-d(t)S 

for all i > To. For any solution x(t) of ([5]) with x(To) = S(Tq) we have by compar- 
ison S(t) < x(t) for all t > T . Let e 2 > satisfy 

£2 4 min {^}- (34) 

Since x*(t) is globally uniformly attractive, there exists T\ > To such that, for all 
t > T\ , we have 

\x(t)-x*{t)\ < e 2 . 
Therefore, for all t > T\, we have 

S(t) <»(t) <x*(t)+£ 2 - (35) 
Thus, by the second equation in (JTJ) and by (Q| we conclude that 
/'(*) _ <p(t,S(t),R{t),Q{t),I(t)) 



I(t) I{t) 



( 7 (t) + a(t) + d(t) + ai (t)) 



< sup !?MMM_ (7(t)+ , (i)+4i)+(ll(t)) 

0<r<A +< 5 1 T 

= &s,« 1 M(t)) 

and finally 

Therefore, using ([331) and J?]) we have, 

lit) = I(T l )e^T 1 b s( T 'S(r))dT < j,j,^ e f^ i b s (T,x*(T)+e 2 )+2e 2 KdT 

By ((321) and dS} we conclude that 

I(t) < /(T 1 )e( Tl -* ) (f- 2 ^^)+ w < I( Tl )e-£ t+N (£ +1 ) 
and using f[33]) and (f34]) we conclude that 

I(t) < HT 1 ) e {^K-N)t+2e.KT 1 < ^ e _£ t +2 E2 im _ 

In any case we have lim /(i) = and this proves that the infectives go to extinc- 

t— >+oc 

tion. 
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Next, still assuming that i? ei i(A) < for some A > or that i? e ,2 < 0, we 
let (S (t) , R(t) , Q (t) , I (t)) be a solution of ((l) with non-negative initial conditions 
and (So(t), 0, Ro(t), Qo{t)) be a disease-free solution of ([T]) with non-negative initial 
conditions. By[m]) in Proposition [1] we can assume that for some S g]0, M[, we have 
(S(t),I(t), R{t),Q(t)), {S (t),0, R (t),Q (t)) G A 4l for t > T . Thus 

S (t),R (t),Q (t), S(t),R(t), Q(t) < A/d + Si, 

for t > Tq. By Theorem [2] we have lim I(t) = 0. Therefore given e > there 

t— >+oo 

exists T e > T such that I(t) < e, S(t),R(t), Q(t) < A/d + Si for t > T e . Thus 

A(i) - d(t)S > S"(t) > A(t) - <f(i)S - if (A/d + Si)e, 

for all t>T e . By comparison we get 

xi(t) > S(t) > x 2 (t) (36) 

where xi(t) is the solution of §5§ and 2:2 (t) is the solution of (|18|) with xi(T e ) = 
x 2 (T £ ) = S(T £ ). ByEvJl in Proposition H we obtain, for all t > Ti, 

K(A/d + 5i)e ,„„. 
Mt)-a:i(t)l<— ^ — , (37) 

Since So(t) is a solution of ([5]), it is globally uniformly attractive, by|mj in Propo- 
sition [2j Thus, there is T2 > T £ such that 

\xi(t) - S (t)\ <e (38) 
for all t > To. Bv (|36 ]) . ([37 ]) and dHJ) we conclude that 

K(A/d + ^Je ^ c m f Jf(A/d + 5i) 
d 

for all t > T2, and also, by (|3"6")l and 

S(t)<Sb(*) + e (40) 

for all t > T2. By (|39|) and (l4Cfl) we obtain, since e > can be made arbitrarily 
small by taking t > T e , we have 

lim \S(t)-S (t)\ = 0. (41) 

Let T 3 > T e be such that \I(t) - I (t)\ < £ for every t > T 3 and write u(t) = 
Q(t) — Qo(t). By the third equation in (JTJ), we get 

-C,u(t) < u'{t) < eoi - (u{t), 

where £ = sup s>0 (<i(s) + Q!2(s) + e(s)) and 01 = sup t>0 a(t). We obtain 

w(T 3 )e- c(t - T3) < u(t) < uiT^e"^-^-^ . 

and thus 

lim \Q(t) -Qo(*)| = lim u(i) = 0. (42) 

Let T4 > Ti be such that — Qo(i)l < e f° r everv t > T4 and write = 

— Ro(t). By the fourth equation in (TTJ) and since Io(t) — and for all t > 0, 
we get 

—eg — diw(t) < w (t) < ep — dw(t), 
where p — sup 4>0 (7(i) + s(t)), q — sup t>0 e(t) and di = sup t>0 d(t). We obtain 

w{T A )e- d ^ T ^ < w(t) < wiTje-W-^-* . 



S(t) > x 2 (t) > xi(t) - y * j L > > S (t) - ^ v 1 d u +l)e, (39) 
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and thus 

lim \R(t) - Ro(t)\ = lim w(t) = 0. (43) 

t— S- + 00 t— > + 0O 

By (l4"Tj) . (1421) . (|43l) and since all the infectives go to extinction, we conclude that 
any disease-free solution (S'o(i), 0, Ro(t), Qo(t)) with nonnegative initial conditions 
is globally attractive. We obtain [51 in Theorem [1] and this concludes our proof. 
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A NONAUTONOMOUS EPIDEMIC MODEL WITH GENERAL 
INCIDENCE AND QUARANTINE 



CESAR M. SILVA 



Abstract. We obtain conditions for eradication and permanence of the in- 
fectives for a nonautonomous SIQR model with time-dependent parameters, 
that are not assumed to be periodic. The incidence is given by functions of all 
compartments and the threshold conditions are given by some numbers that 
play the role of the basic reproduction number. We obtain simple threshold 
conditions in the autonomous, asymptotically autonomous and periodic set- 
tings and show that our thresholds coincide with the ones already established. 
Additionally, we obtain threshold conditions for the general nonautonomous 
model with mass-action, standard and quarantine-adjusted incidence. 



1. Introduction 

Quarantine is one of the possible measures for reducing the transmission of dis- 
eases and has been used over the centuries to reduce the transmission of human 
diseases such as leprosy, plague, typhus, cholera, yellow fever and tuberculosis. 

The concept of quarantine dates back to the fourteenth century and is related 
to plague. In fact, in 1377, the Rector of the seaport of the old town of Ragusa 
(modern Dubrovnik) officially issued a thirty-day isolation period: ships coming 
from infected or suspected to be infected sites were to stay at anchor for thirty 
days before docking. A similar measure was adopted for land travelers but the 
period of time was enlarged to forty days. It is from the Italian word for forty - 
quaranta - that comes the term quarantine [9] . 

Isolation is still nowadays an effective way of controlling infectious diseases. For 
some milder infectious diseases, quarantined individuals could be individuals that 
choose to stay home from school or work while for more severe infectious diseases 
people may be forced into quarantine. 

To study the effect of quarantine in the spread of a disease, we consider a SIQR 
model. Thus, we divide the population into four compartments: the S compartment 
corresponding to not infected individuals that arc susceptible to the disease, the I 
compartment corresponding to individuals that are infected and not isolated, the Q 
compartment corresponding to isolated individuals (Q stands for quarantine) and 
the R compartment including the recovered and immune individuals. We assume 
in this paper that the quarantined individuals are perfectly isolated from the others 
so that they do not infect the susceptibles and that the infection confers permanent 
immunity upon recovery. 
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A major concern when studying epidemic models is to understand the asymptotic 
behavior of the compartments considered, particularly the infective compartment. 
Frequently, in autonomous epidemic models, one of two situations occur: there is a 
disease-free equilibrium that is the unique equilibrium of the system or the disease- 
free equilibrium coexists with an endemic equilibrium. The basic reproduction 
number is a fundamental tool to obtain the behavior for a given set of parameters: 
if the basic reproductive number is less than one the disease-free equilibrium is 
locally asymptotically stable and if this number exceeds one the disease remains in 
the population. 

Several authors studied autonomous epidemic models with a quarantine class. 
Namely, Feng and Thieme [8] proposed a SIQR model for childhood diseases and 
showed that isolation can be responsible for the existence of self-sustained oscilla- 
tions. Their model includes a modified incidence function - the quarantine-adjusted 
incidence - where the number of per capita contacts between susceptibles and infec- 
tives is divided by the number of non isolated individuals. Some dynamical aspects 
of this SIQR model were studied by Wu and Feng [19]. Namely, these authors 
studied this model through unfolding analysis of a normal form derived from the 
model and found Hopf and homoclinic bifurcations associated to this unfolding. A 
generalized version of Feng and Thieme's model, allowing infected individuals to 
recover without passing through quarantine and also disease related deaths, was 
proposed by Hethcote, Zhien and Shengbing in [5], where this model was compared 
to SIQR models with different incidences. A similar analysis is also undertaken in 
that paper for SIQS models with different incidences. The authors concluded in 
that work that the SIQR model with quarantine-adjusted incidence seems consis- 
tent with the sustained oscillations that are observed in real disease incidences but 
do not match every aspect of the observed data. 

Other types of epidemic models with quarantine were studied in the literature. 
A SEIQR model was considered by Gerberry and Milner in [4] where the authors 
studied the dynamics and used historical data to discuss the adequacy of the model 
in the context of childhood diseases. The global dynamic behavior of a different 
SEIQS epidemic model, with a nonlinear incidence of the form f3SI, was analyzed 
by Yi, Zhao and Zhang in |16j . In another direction, Arino, Jordan and van den 
Driessche pQ studied the effect of quarantine in a general model of a disease that 
can be transmitted between different species and multiple patches (though in their 
setting quarantine is present only in the form of travel restriction between patches). 
A model with twelve different classes was proposed by Safi and Gumel [T5] to study 
the impact of quarantine of latent cases, isolation of symptomatic cases and of 
an imperfect vaccine in the control of the spread of an infectious disease in the 
population. The same authors studied a delayed model including a quarantined 
class [14]. 

We emphasize that all the above models are autonomous models. However, it is 
well-known that fluctuations are very common in disease transmission. For instance, 
the opening and closing of schools tend to make the contact rates vary seasonally in 
childhood diseases [3] and periodic changes in birth rates of the population are also 
frequent. It is therefore natural to consider time-dependent parameters in epidemic 
models and define threshold conditions for extinction and persistence of the disease. 

Our objective in this paper is to consider a family nonautonomous SIQR models 
with a large family of contact rates given by general functions that are, in general, 
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time-dependent and far from bilinear, and obtain numbers that play the role of 
the basic reproductive number in this setting. Related results for nonautonomous 
SIRVS models were obtained in JT2J [2U] , thought the incidences considered in those 
papers are linear in the infectives while here we assume that the incidences are, in 
general, nonlinear functions of all compartments. 

Note that in p~8l [13] some methods were developed to obtain thresholds for 
general periodic models. These methods apply to our setting when we assume the 
coefficients to be periodic. We emphasise that this is not necessarily our case, since 
our coefficients need not be periodic. Additionally, our thresholds are given by 
explicit formulas (see section [3]). 

The structure of the paper is as follows: in section[2]we introduce the model and 
the main definitions; in section [3] we state our main Theorem; in the next sections 
we look at some particular cases of our model, namely, in section|4]we particularize 
for the autonomous model, in section [5] we see that the asymptotically autonomous 
model has the same thresholds as the limiting autonomous model, in section [6] we 
obtain thresholds for the periodic model with constant death and recruitment rates 
and, in section [7] we consider the general nonautonomous model with mass-action, 
standard and quarantine- adjusted incidences; in section [8] we prove some auxiliary 
results and the main theorem; finally, in section [5] we present some simulations that 
illustrate our result. 



2. Generalized nonautonomous SIQR model 

We consider the following nonautonomous and general incidence SIQR model 
that generalizes the three autonomous SIQR models in [S], 

' S' = A(i) - ip(t, S, R, Q, I) - d(t)S 
V = <p(t, S,R,Q,I)- [ 7 (t) + <7(t) + d(t) + at (t)] I 
Q' = a(t)I - [d(t) + a 2 (t) + e{t)}Q ' U 

_R' = j(t)I + e{t)Q -d{t)R 

where S, R, Q and / correspond respectively to the susceptible, recovered, quar- 
antine and infective compartments; A(t) is the recruitment of susceptible (births 
and immigration); (p(t, S, R,Q, I) is the incidence (into the infective class); d(t) is 
the per capita natural mortality rate; j(t) is the recovery rate; a(t) is the removal 
rate from the infectives; a\(t) is the disease related death in the infectives; 0:2 (t) is 
the disease related death in the quarantine class; and e(t) is the removal rate from 
the quarantine class. We will assume that (p, A, d, 7, a, ot\, a 2 and e are continuous 
bounded and nonnegative functions on and that there are u>d,UA > such that 

d~ d >0 and AJ A >0, (2) 

where we are using the notation 

h~ =liminf / h(s) ds and hz = limsup / h(s)ds, 

h J t h t^+oo J t 

that we will keep on using throughout the paper. For bounded h we will also use 
the notation 

hs=swph(t) and hi = inf (t). 
Next we state some simple facts about problem (JTJ) . 
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Proposition 1. We have the following: 

i) all solutions (S(t), I(t), Q(t), R(t)) of JTJ with S(t ) > 0, I(t ) ^ 0, Q(t ) > 
and R(t )) ^ verify S(t) ^ 0, I(t) ^ 0, Q(t) ^ and R(t)) ^ for all t ^ t ; 

ii) all solutions (S(t), I(t),Q(t), R(t)) of JTJ with S(t ) > 0, J(t ) > 0, Q(t ) > 
and R{t )) > veri^ S(t) > 0, I(t) > 0, Q(t) > and R(t)) > for all t ^ t ; 

Hi) If (S(t), I(t),Q(t), R(t)) is a solution of JTJ with nonnegative initial conditions 
then there is a constant K > such that 

limsup5(t) + I{t) + Q{t) + R(t) sc K. 

t—t+oo 

iv) There are constants C,T > such that, if I(t) < (3 for all t > T then 

Q{t) < Cf3 and R(t) < C[3 (3) 
for all t sufficiently large. 

Proof. Properties [jj , [nj are a simple consequence of the nonnegativeness and the 
boundedness of the parameter functions. Using J3J, property [mj follows easily. 
Assuming that I(t) < (3 for all t > T, by the third equation in JTJ, we get 

Q'{t) < (3a{t) - [d(t) + a 2 {t) + e(t)]Q(t) 
and thus, by J3J, for some T > T sufficiently large and Co = ^f, we have 

Q{t) = e - & d ( s )+ a 2(s)+e(s) ds Q( T j + e -fl d(s)+ a2 { s )+e{s) ds ^ 

JT 

^ e -C (t-T-uj d ) q(t) + Pas [ e - c °^- u - w ^ du 

JT 

^ e -Co(t-T-u d) R + ffosf^ A _ e -C 0{ t-T)\ 

Cq V / 

and thus, taking C x = cs e Co " d /(2C ), we have Q(t) < C x f3 for t sufficiently 
large. A similar argument applied to the third equation in JTJ shows that there is 
a constant Ci > (independent of (3) such that R(t) ^ 6*2/3 for t sufficiently large. 
Choosing C — max{Ci, C2} we obtain property HvJ . □ 

Consider the auxiliary equation 

x' = A(i) - d(t)x, (4) 

Proposition 2. We have the following: 

i) Given to 0, all solutions x(t) of equation (fTJ with initial condition x(to) 

are nonnegative for all t ^ 0; 
ii) Given to ;J 0, aZ/ solutions x(t) of equation J4j ure'i/i initial condition x(to) > 

are positive for all t ^ 0; 
mj Given to ^ and a solution x(t) of equation J3J wi£/i initial condition x(to) > 0, 

£/iere are M > and t\ smc/i £/iai :c(i) ^ M /or all t^ti; 
iv) Each fixed solution x(t) of (TTJ with initial condition x(fo) ^ is bounded and 

globally uniformly attractive on [0, +oo[; 
v) There is a constant D > such that, if x(t) is a solution of (T4J and x(i) is a 

solution of the system 

x' = A(t) - d(t)x + f(t) (5) 
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with x(to) = x(to) then 

sup |x(£) - x(t)\ < D sup \f(t)\. 

Proof. For to 0, the solution of ((4]) with initial condition x(to) — xo > is 
given by 

x(t) = c~ d(s) ds x + f e~ /- d < s > ds A(u) du. (6) 



Since £ i— > A(£) is nonnegative, equation ([5]) immediately implies [I]) and|n]). 
Since t h-> is bounded, by and ([6]) and ([2]) we obtain 

x(t) = e~ f*-»A d{s) ds x (t - lja) + / e - £ ds A(«) 



t-Ul A 



>e~ dsUA / A(u)dw 

•/ t — UJA. 



for all £ sufficiently large. Thus we have lm|) . 

By ©, since £ M> A(£) is bounded, we get for all £ sufficiently large 

x(t) = e~ & d{s) ds xo + f e~ £ d(s) ds A(u) 



< x + A s / exp ( -- — - - 1 J 4r rfu 



'to 



w d / 2 



<* + ^V-/ 2 A s , 



and we conclude that each solution x(£) of ((4]) with initial condition x(to) 
is bounded. Additionally, if x, x\ are two solutions of (j4|) with initial condition 
xo = x(£o) ^ and xi.o = xi(£o) we have 

|z(t)-zi(t)|<e-£ dW *|so-zi,o|, 

and we get |Iv]) . 

Finally, if x(£) is a solution of (0} and x(£) is a solution of (J5J) with x(£o) = 
x(£o) then, setting u(t) = x(t) ~ x(£), we obtain the problem 

v! = -d{t)u + /(£), 

with u(£o) = 0. The above problem has solution 

u{t)= / c-ti d ^ ds f{u)du 



to 



ft 

< sup|/(£)| / e -X dis)ds du 

t>to Jt 



<£>sup|/(£)| 

t>t 



where D = ^ exp f ^-J . Thus, we obtain |Vj|. □ 
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For each 8,5 £ Rg with 5 > 8 define the set 

A e ,«5 = {(xi,x 2 ,x 3 ,x i ) el 4 : < xi < 5 A < Xi < 8, i = 2,3,4}. 

We note that every solution (S(t),I(t), Q(t),R(t)) of our system stays in the region 
Ao,k (where if is givel by[mj in Proposition [IJ for every t 6 R^ sufficiently large. 

We need some additional assumptions about our system. Assume that: 
Hf ) the functions x i-> x, y, it), z) are non decreasing and ip(t, 0, j/, u>, z) = 0; 
H2) given 8 > there is AT# > such that 

\ip(t,xi,y,w,z) - ip(t,X2,y,w,z)\ s$ Ar |xi - x 2 \z, 

for each (t, Xi, y, w, z), (t,x 2 ,y,w,z) Gljx Ab,k\ 

. - - In ,sl .,,,,1 ; - 1 ( 



H3) There is N > such that, for every 5 € ]0, K], y,W,z € ]0, 5] and t e lj we 



have 



(p(t,x,0,0,T) < ip(t,x,y,w,z) <N 



0<r<35 T 

H4) For every y,w, z G ]0, A"] we have 

<p(t,x,y,w,z) <p(t, x, 0, 0, 5) 
^ limsup < +oo. 

Note that, for contact rates that are not dependent on y and u>, is always 
verified. Note also that, if for each 8 s]0, AT] there is AT# > such that 

dip 

— (t,x,y,w,z) ^ AT e z, 

for all (t, x,y,w, z) £ Rj x A^e, then holds. Now, we need to make some 
definitions. We say that 

i) the infectives I are permanent in system ([T]) if there exist constants < m\ < 
m 2 such that 

mi < liminf I(t) ^ lim sup /(i) < m 2 , 

t~¥O0 f-J-OO 

for every solution (£(*), I(t), R(t), V(t)) of ([I]) with S (t), I {t), R {t),V (t) > 
(note that we require that mi and m 2 are independent of the given solution 
with positive initial condition in /); 

ii) the infectives / go to extinction in system (fTJ) if lim I(t) = for all solutions 

t— >oo 

Of ©. 

3. Statment of the results 
For each solution x(t) of (U]) with x(0) > 0, define the function 

bg ^ x{t)) = <p(t,x(t),0,0,6) _ (7(t) + a(t) + m + Qi(t)) 



and the numbers 



r p (X) = liminf / liminf bs (s, x(s)) ds, 

*->+oo J. 8— >0 
i-t+X 

r e (A) = limsup / lim sup b$(s, x(s)) ds, 

t— > +oo Jt S— yO 

R p (\)=e r * w and R e (X) = e r = (A) . 
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Contrarily to what one could expect, the next technical lemma that will be proved 
in Section |8] shows that the numbers above do not depend on the particular solution 
x(t) of (U} with a;(0) > 0. The numbers R P (X), R e (X) play the role of the basic 
reproduction number (defined for autonomous systems and also called quarantine 
number) in the non-autonomous setting. 

We have the following result whose proof will be given in section [8] 

Lemma 1. We have the following: 

1. Let e > be sufficiently small and < 9 ^ K . If 

a, b€]9,K[ and a — b<e, 

then 

\b s (t,a)-b 5 (t,b)\<sK e . (7) 

2. The numbers Rp{X) and R e (X) are independent of the particular solution 
x(t) with x(Q) > of fll}- 

Note that, by ([7]) we have 



lim inf b$ (t, a) 

5^0 



lim inf b$ (t, b) 

5^0 



and 



lim sup bg (t , a) — lim sup bs (t, b) 
<5->o <5^o 



< eK e 



< eK e . 



(8) 



(9) 



We now state our main theorem on the permanence and extinction of the infec- 
tives in system (fTJ). A proof of this result will be given in section [8] 

Theorem 1. We have the following for the system (fT|). 

1. If there is a constant A > such that R P (X) > 1 then the infectives I are 
permanent. 

2. If there is a constant A > such that i? e (A) < 1 then the infectives I 
go to extinction and any disease free solution (S(i),0, R(t),Q(t)) is globally 
attractive. 



4. Autonomous model 



In this section we are going to consider the autonomous setting. Namely, we are 
going to assume in system ([l} that A,d, 7, a, a.\,E and c*2 are constant functions, 
that A, d > 0, and that (p is independent of t. We obtain the autonomous system 



'S' = A-<p (S,R,Q,I) -dS 

r = <p {s,R,QJ)-h+^ + d- 

}' = o-I - [d + a 2 +e]Q 
^R' = 7/ + eQ - dR 



(10) 



where tpo : (H,}) 4 — > E, is continuous, nonnegative and satisfies t{T]), an d EE])- 
In this setting we have that the auxiliary equation ((3]) admits the constant solution 
x(t) = A/d. Thus we obtain 

<A,(AM0,0,<5) 



lim inf ■ 

5^0 



(7 + a + d + ai) 



A, R P (X) = c 
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and 



lim sup (7 + (T + a + ai) 

5^0 o 



A, i? e ,i(A)=C 



It is now easy to establish a result that is a version of Theorem [T] in the particular 
case of autonomous systems. Define 



<5^o 5(7 + cr + d + ai) 5_> 5(7 + cr + (i + Q;i)' 

For any A > 0, we have i? p (A) > 1 if and only if Rp UT > 1 and R e (X) < 1 if and 
only if R^ UT < 1 and this implies the following result. 

Theorem 2. We have the following for the autonomous system (|10|) . 

1. If Rp UT > 1 then the infectives I are permanent; 

2. If R^ UT < 1 then the infectives I go to extinction and any disease free 
solution is globally attractive. 

Note that, if the incidence has the particular form 

<Po(S, R, Q, I) = ip(S, R, Q)g(I)I, (11) 

where g is a continuous, bounded and nonnegative function, then 

ip(A/d, 0, 0) liminf g(S) ^(A/d, 0, 0) limsup3(<5) 

R AUT = S±0 ^ R AUT_ 6^0 



p 7 + a + d + ol\ ' 7 + a + d + a.\ 

and we conclude that the infectives are permanent (resp. go to extinction) if 
, . 7 + a + d + «i 7 + o- + d + a x 

W ^ > m/d,o,o) (resp - hl L s up5(5) < #Mo,o) } ' 

Naturally, if 

7 + cr + d + ai 



G 



liminf g(S), limsupg(<5) 



^(A/d,0,0) 

we get no information about the asymptotic behavior of the infectives. 

Assuming now that g is constant, we have Rp UT — R^ UT and we can recover 
the results obtained in J5]. In fact, the autonomous SIQR model with mass-action 
incidence (<po(S, R, Q, I) = /3SI), the autonomous SIQR model with standard in- 
cidence (ipo(S,R,Q,I) — (3SI/{S + I + Q)) and the autonomous SIQR model 
with quarantine-adjusted incidence (<po(S, R,Q, I) = (3SI/(S + I + Q + R)) are 
all in the conditions we have assumed and our threshold values for those model 
are Rp UT — Rf UT — /3A/[d(7 + a + d + ai)] for mass-action incidence and 
rAut _ rAut _ _|_ a _)_ d _|_ ai ) f or standard and quarantine-adjusted inci- 
dences. For these models, the numbers Rp UT — Rf UT coincide with the quarantine 
reproduction number obtained in [5]. Note that, to put the referred autonomous 
models considered in [5] in our context, we need to have the standard incidence and 
the quarantine-adjusted incidence continuous in A^o- This constitutes no problem 
since the functions tpo associated to these incidences can be extended to a con- 
tinuous function in A^Oj setting tpo(0, R, 0, 0) = for the standard incidence and 
t/?o(0, 0, 0, 0) = for the quarantine-adjusted incidence. 

Set ip(S, Q,R) = S and g(I) = withp, q > 0, a ^ 0, in dHJ). This family of 

contact rates was considered for instance in [10 , 6 . Our threshold conditions show 
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that if p < 1 then the disease is permanent (independently of the other parameters) 
and if p — 1 , setting 

V»(A/d,0,0) 

IX = : , 

7 + <j + a + a.\ 

the disease is permanent if R > 1 and goes to extinction if i? < 1. 

5. Asymptotically autonomous model 

In this section we are going to consider the asymptotically autonomous SIQR 
model. That is, additionally to the assumptions on Theorem [TJ we are going to 
assume for system ([1]) that there is a continuous function (po such that 

lim <p(t, S, R, Q, I) = (f (S, R, Q, I) 

t— >oo 

for each (S, R, Q, I) £ R 4 and that the time-dependent parameters are asymptot- 
ically constant: j3(t) -> f3, A(t) -4 A, d(t) -> d, 7(f) -4 7, cr(t) -> ex, ai(f) -> a\, 
oc%(t) — > oli and e(t) — >■ £ as t — > +00. Denoting by F(i, S 1 , i?, Q, /) the right hand 
side of ([1]) and by Fq(S, R, Q, I) the right hand side of the limiting system, that is 
of (fTT)|) , we also need to assume that 

lim F(t,S,R,Q,I)=F (S,R,Q,I), 

with uniform convergence on every compact set of (R^) 4 , and that (S, R, Q, I) 1— > 
F(t, S, R, Q, I) and (S, R, Q, I) M- F (S, R, Q, I) are locally Lipschitz functions. 

There is a general setting that will allow us to study this case. Let / : Rxi?" — »■ R 
and /o : R" -4 1 be continuous and locally Lipschitz in R™. Assume also that the 
nonautonomous system 

x' = f(t,x) (12) 
is asymptotically autonomous with limit equation 

x 1 = f (x), (13) 

that is, assume that f(t, x) — > fo(x) as t — > +00 with uniform convergence in every 
compact set of R™. The following theorem is a particular case of a result established 
in [IT] (for related results and applications see for example [H [7]). 

Theorem 3. Let <I>(t, to, xq) and <p(t,to,yo) be solutions of (|12p and (I13[) respec- 
tively. Suppose that e £ R™ is a locally stable equilibrium point of (|13p with attrac- 
tive region 

W(e) = \yeR n : lim <p.(t,t ,y) = e\ 

and t/iat W$ PI W(e) ^ 0, where W$ denotes the omega limit of <I>(t, to, xo). TTien 
lim <I>(t, to, £0) = e. 

t— >+oo 

Since (R + ) 4 is the attractive region for any solution of system (jTU)) with initial 
condition in (R + ) and the omega limit of every orbit of the asymptotically au- 
tonomous system with /(to) > is contained in (R + ) 4 , we can use Theorem [5] and 
Theorem [3] to obtain the following result. 

Theorem 4. We have the following for the asymptotically autonomous systems 
above. 

1. If Rp UT > 1 then the infectives I are permanent; 
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2. If Rf UT < 1 then the infectives I go to extinction and any disease free 
solution is globally attractive. 

6. Periodic model with constant natural death and recruitment 

In this section we are going to consider a periodic SIQR model. Additionally 
to our assumptions on the function tp and the parameter functions in Theorem [TJ 
we are going to assume in system ([1]) that A(i) = A and d(t) = d are constant 
functions, that there is a T > such that ip(t, S, R, Q, I) — ip(t + T, S, R, Q, I) and 
that the remaining time-dependent parameter functions are periodic functions with 
period T. We have in this case the constant solution x{t) = A/d and therefore 

r p (A)=liminf / + lim inf ^ °' °' ^ - (7(5) + <r(s) + d + a 1 (sj) ds 



t->+oo L <5->0 



r liminf f(s,\/d o A s) ds _ +d 

Jo s^o S 



and similarly 



r e (A) = / hmsup as — (7 + a + a + (Xi)T, 

Jo s->o 

where / denotes the average of / in the interval [0,T]: / = ^ J Q T f(s) ds. Define 



I [ T , , ip(s,A/d,0,0,6) 

and 



R P P ER W = ^ Hniinf "'V s 

p w T J Q s^o 5(j + a + d + ai ) 



K ER W = U hmsup M^fc 



T Jo T^cT T(7 + CT + rf + ai) 

Theorem 5. We have the following for the periodic system with constant recruit- 
ment and death rates. 

1. If Rp ER > 1 then the infectives I are permanent; 

2. If R P < 1 then the infectives I go to extinction and any disease free 
solution is globally attractive. 

If we assume in the periodic model above that the incidence has the particular 
form ip(t, S, R, Q, I) = /3(t)SI with j3(t + T) = /3(t), then we obtain 

R PER (\) = R PER {\) = -. ^ r (14) 

and if we assume in that model that ip(t, S, R, Q, I) = f3(t)SI/(S + R + Q + I) or 
that <p(t, S, R, Q, I) = /3(t)SI/(S + R + I) with (3(t + T) = (3(t) we get 

R PER (X) - R PER W = — j^— ■ (15) 
7 + a + d + ai 

If the parameter functions are all constant, by (TUT) and (IT51) we obtain again the 
thresholds for the autonomous system with mass-action, standard and quarantine- 
adjusted incidence. 

The thresholds (fT4")) and (fT5")) can also be obtained using the methods developed 
in jTS] [T3] for general periodic epidemic models, that constitute periodic versions 
of the general autonomous model considered in [T7] . 
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7. General model with mass-action, standard and 
quarantine-adjusted incidence 

In this section we are going to consider a nonautonomous SIQR model with 
mass-action, standard and quarantine-adjusted incidence. Additionally to our as- 
sumptions on the function ip and the parameter functions in Theorem [TJ we are 
going to consider the particular cases of mass-action, standard and quarantine- 
adjusted incidence. We will also assume that natural death and recruitment are 
constant. We have in this case the constant solution x(t) = A/d. 

For mass-action incidence, tp(t, S, I, Q, R) = f3(t)SI, we have 



Una inf 

t— >+oo 

A 



- il{s)+a(s) +d + a 1 (s))ds 



^ A 



^x ~hx + ct a + d + (ai) A 



and similarly 



Define 



r e (A) > A [7+ + a+ + d + ( ai )+]) . 



i? 



SIM 



(A) 



A/3, 



and 



7? 



SIM 



(A) 



d [7 A + °a + ("i)a1 
A/?a + 



Theorem 6. FFe have the following for the general non- autonomous system with 
mass-action incidence and constant recruitment and death rates. 

1. If there is A > such that Rp IM (X) > 1 then the infectives I are permanent; 

2. If there is A > such that i?f /M (A) < 1 then the infectives I go to extinction 
and any disease free solution is globally attractive. 

For standard and quarantine-adjusted incidence, respectively 

[3(t)SI j _ u a T „ m (3{t)SI 



<p(t,S,I,Q,R) = 



S+I+R+Q 



and ip(t,S,I,Q,R) 



S + I 



we have 
r P (A) 



lim inf 



* +A A/d 
P(t) lim inf- 7 



s->o A/d + 5 

<M^A ~[7a + ct a +^+(«i)aO 
and similarly 

r e (A) > A(/3+-[ 7 + + ( T++d+( ai )+]) 



(7(5) + cr(s) + d + ai(s))ds 



Define 



and 



Rp^ A (X) = 



R s p / QA (X) = 



7a 



d+ (ai); 



It 



7 + + ( j++d+(a 1 )+ 
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Theorem 7. We have the following for the general non- autonomous system with 
mass-action incidence and constant recruitment and death rates. 

1. If there is A > such that R P ^ A (X) > 1 then the infectives I are permanent; 

2. If there is A > such that Re^ A (X) < 1 then the infectives I go to extinction 
and any disease free solution is globally attractive. 

8. Proof of Theorem Q] 

8.1. Proof of Lemma [TJ Assume that e > 0, that < 6 < K , that a, b e]9,S[ 
and that a — b < e. We have, by E(2]), for every < r < S, 

\<p(t, a, 0, 0, 5) ~ tp(t, b, 0, 0,(5)1 < K e \a - b\8 < Kgbe. 

Therefore, 

<p(t,a,0,0,S) v „ <p(t,b,0,0,S) „ y(M,0,0,3) , „ 

K g e < < 1- K e e, (16) 

and, adding and subtracting j(t) + a(t) + d(t) + ai(t), we get ©. 

We will now show that in fact r p (X), r e (X) are independent of the particular 
solution x(t) of ((4} with x(0) > 0. In fact, by Ov} in Proposition [2j for every 
e > and every solution Xi(t) of (0| with xi(0) > there is a T e > such that 
\x(t) — xi(t)\ < e for every t T £ . Choose > such that, for e < Ek and 
t ^ T s , we have x(t), x\ (£), x(t) ±ee]S, if [, for some > that depends on x and 
Xi. Lcting a = x(t), b = x\(t) we obtain by (jSJ 

/■i+A /-t+A 

/ lim inf bg (s, £i(s)) c£s — XKge < / lim inf bs(s, x(s)) ds 
Jt Jt 

r t+x 



< / lim inf bg(s, x\ (s)) ds + XKge 
Jt 

every < e < ek, 

j-Z+A pt + X 

liminf / liminf bg(s, x(s)) ds — liminf / liminf bg(s, Xx(s)) ds 

t— >+oo J. 8— >0 t— !-+oo J t 5—^0 



for every t ^ T e . We conclude that, for every < e < Sk, 

ft+X rt+X 



< XK g e, 



it 

and thus 

/•t+A pt+X 

liminf / liminf bs(s, x(s)) ds = liminf / liminf bg(s, Xi(s)) ds. 

t—t+oo J t S— s-0 t— >+oo J t 

Letting 5 — > we conclude that r p (X) (and thus R P (X)) is independent of the 
chosen solution. Using (|9]), the same reasoning shows that r e (X) (and thus R e (X)) 
are also independent of the particular solution. This proves the lemma. 

8.2. Proof of Theorem [TJ Assume that R P (X) > 1 (and thus r p (X) > 0) and let 
(S(t)J(t),R{t),Q(t)) be any solution of © with S(T ) > 0, I(T ) > 0, R(T ) > 
and Q(Tq) > for some To ^ 0. By[m]) in Proposition [TJ we may assume that for 
t Ss T we have (S(t), I(t), R(t),Q(t)) e A , K - 

Since r p (X) > 0, using (jHJ we conclude that there are constants 5i s]0, K[ and 
> such that 

r t+X 

inf 6 5 (s,x*(s) -e)ds > 6 (17) 

t 0<<5<3<5i 

for all £ ^ sufficiently large, say £ Ti, and e > sufficiently small, say for 
< £ < e, where x* is any fixed solution of (|4]) with x*(0) > 0. Note that, by[u| in 
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Proposition [2] we have x*(Tq) > and by by|m| in Proposition [2] we may assume 
(eventually using a bigger Ti ) that x* it) ^ 9 for some 6 > and i ^ Ti. 
Define 

£o = min {4^A' £ "} and £l=min { 5l '§'2^v}' (18) 

where AT^ is given by t(2]), Af is given by tJ3]), C is given byEy} in Proposition [1] 
and D is given by|vj) in Proposition [2J 
We will show that 

lim sup J(i) > e x . (19) 

t— >+oc 

Assume by contradiction that (|T9l is not true. Then there exists T 2 ^ satisfying 
J(t) < £i for all t T 2 . Consider the auxiliary equation 

x' = A(t) - d(t)x - Ne u (20) 

where N is given by tj3|). 

Let x(t) be the solution of (j4j) with x(Tq) = 5 (To) and let x(t) be the solution 
of J2(| with x(T ) = 5(T ). By (JTSj) and[vj) in Proposition [2] we obtain 

\x(t)-x(t)\^DN El ^^, (21) 

for all t > T . 

According to Gv} in Proposition [21 x* (i) is globally uniformly attractive on Rq . 
Therefore, there exists T3 > such that, for all t ^ T3, we have 

|*(t) -**(*)! <y. (22) 

By t(3l) we have 

sc <p(t, S(t),R(t), Q(t),I(t)) < AT/(i) < A^, 

for all i ^ max{Tb, T 2 }, and thus, by the first equation in ([T]), 

A(i) - d(t)5 > 5' > A(i) - d(t)S - AT £l . (23) 

for all t ^ max{T ,r 2 }. Comparing (J20j and ([23} we conclude that 5(t) > x(f) for 
all f > max{r ,r 2 }. Write T 4 = max{T , Ti, T 2 , T 3 }. By ([21]) and ([221) we have, 
for all t > Ti, 

X*(t) + £ > X*(t) + y ^ 5(t) > 5(4) > x(t) > x(t) - y > X*(t) - £ , 

and thus 

\S(t) - x*(t) + £ | < \S(t) - x*(t)\ + £ < 2£ , (24) 
for all t ^ Ti. By ([24]) and in Lemma Q] we get, for all t > T 4 , 

b 5 (t,S(t)) >b s (t,x*(t)-e )'2Kge . (25) 

For i > sufficiently large, say t ^ T 5 , we have by © and ([IS]) . Q(i) < Cei < Si 
and < Cei < S± and thus I(t),Q(t), R(t) £ ]0,<5i[. We may assume that 

T 5 ^ T 4 . 

According to (T5]) we have £0 < £ and by (fTTj) we get 

inf b s (T,x*(s)-s )ds> (ht-T s )-l)e. (26) 

ji 0<5<35i y A y 
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By the second equation in ([T]) we obtain 
-cp(t,S(t),R(t),Q(t),I(t)) 



I'{t) = 



7(f) 



- ( 7 (t) + a(t) + d(t) + ai(t)) 



J(i) 



and thus, integrating from T5 to t, recalling that I (t) , Q (f) , R(t) <E ]0, <5i [ for all 
t > T 5 and using r£J]), (JUJ, ([26j and (HU, we have 

/(<) = I(T 5 ) e-^s ^(^(s),fl( S ),Q(s) J( S ))//(sM 7 (s)+a(s)+d( S )+ai(s)) 
f* inf fc a (s,S(s))ds 

^/(T 5 )e 50<5<3ai 

>/(T 5 )e^»<'< 3 ^ i(S ' X * (s) " 0) ^ 2K9E0(t " T5) 

> J(T B )e<*- 

> J(T 5 )e* 

and we conclude that 7(t) — s> +00. This contradicts the assumption that J(i) < £1 
for all t ^ T2. From this we conclude that (TT91 holds. 

Next we will prove that for some constant I > we have in fact 



liminf m) > I 

t->+oo 



(27) 



for every solution (S(t),I(t),R{t),Q(t)) with S(T ) > 0, J(T ) > 0, i?(T ) > and 
Q(T ) > 0. By CE7]), for all ^ > A, * ^ T x and < < e, we have 

inf bg(s, x*(s) - 9) ds > N. (28) 

We proceed by contradiction. Assume that (f27j) does not hold. Then there 
exists a sequence of initial values (x n ) n ^, with x„ = (S„, R n , Q n , In) with £„ > 0, 
R n > 0, Q n > and J„ > such that 

liminf I(t,x n ) < — k 

where I(t,x n ) denotes the solution of (UJ with initial conditions S(Ti) = S n , 
J(Ti) = I n , R(Tx) = R n and Q{T X ) = Q n . By $TS]|, given n 6 I, there are 
two sequences (i n ,fc)fcew and (sn,fe)fee»r with 

T\ < S„ t l < t„ : i < S n ,2 < t n ,2 < ■ ■ ■ < S n> k < t n ,k < ■ ■ ■ 



and lim s n ,k = +00, such that 

k— >+oo ' 



and 



n n z 

£ 1 , r/, \ £l 



-~ < < — , for all t e.]s n , k , t n J. 

By the second equation in ([TJ we have 

/'(*, x n ) = I(t, x n )[(ip(t, S(t, x n ), R(t, x n ), Q(t, x n ),I(t, x n ))/I(t, x n )- 
-( 7 (t)+5(t)+d(t) + a 1 (t))} 
^ -( 7 (i) + a(t) + d(t) + ai(t))I(t, x n ) 
> -{is +&s + ds + (ai)s)I(t, x n ), 



(29) 
(30) 
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Therefore we obtain 

y, ■ ? \ dr ^ -(7s + a s + d s + (ai)s)(ifc, n - Sfc,«) 

and thus I(t Kni x n ) ^ I(s fcin ,2; n )e~^ s+CTS+ ' is+ ( Ql ) s )( tfc ."- Sfe -). By (JUJ) we get 
I ^ e -(7s+<Ts+ds+(ai)s)(U, n -sfc, n ) an d therefore we have 

log n 

tk,n - S k n ^ ; — — > + OO (31) 

7s + ct s + d s + (ai)s 

asn-> +00. 

According to (|3"0|) . for all < G]sfc.n, ifc,n[, we have I(t,x n ) < £\jn < E\ and we 
conclude that ip(t,S,R,Q 7 I) < Nei for all t ejs^n, Thus by (TTJ) we obtain 

S' ^ A(t) - d(t)S - Nei 

for all t e]s k ,n,tk,n[- 

By comparison we have S(t,x n ) ^ x{t) for all t S]sfc jn , ife, n [, where cc is the 
solution of (f2T)l) with x(sk, n ) = S(sk,n> x n)- Using (p?Tj) we obtain, for all t G 

, (32) 

where x(t) is a solution of (|4]) with x(s/s )n ) = S(sk, n , %n)- Since x*(t) is globally 
uniformly attractive, there exists T* > 0, independent of n and k, such that 

|5(t)-.**(t)|^|, (33) 

for all t ^ Sk, n + T*. By (|3T|) we can choose B > such that i„ jfc - s nik > \ + T* 
for all n ^ £>. Given n ^ B, by (|28|) . (|32|) and (|33l) and by the second equation 
in ([T]) we get 

% = /(**,»,*») > '(•*,« +T*,x„)e^>* li ^o fb5(r ' S(r ^ ))<ir 

£l / . fc '", T ,» lim inf 67. 4 . (r,£c* (r)— e ) dr £\ 

^ — e Sfc ." +J s ->-° > — . 

This leads to a contradiction and establishes that lim inf I(t) > £ > 0. Thus, [I] in 

t— >+oo 

Theorem Q] is established. 

Assume now that R e (X) < 1 (and thus r e (A) < 0) for some A > 0. Then we can 
choose 9 > and To > such that 

t+A 

b s (s,x*{s) +e)ds < -9, (34) 

for all t ^ To and e > sufficiently small, say e < e. By [m} in Proposition Q] we 
can assume that (5(t), 7(t), R(t),Q(t)) E A ,k for £ > T . 
We will show that 

limsupJ(i) = 0. (35) 

By (fTJ) and f{|]) we have 

S'(t) ^ A(t) - d(t)S 
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for all t To. For any solution x(t) of ((4]) with x(Tq) = S(Tq) we have by compar- 
ison S(t) ^ x(t) for all t ^ To. Since x*(t) is globally uniformly attractive, there 
exists £\ < e and T\ ^ To such that, for all t > T-y, we have 

|x(i) - X*(t)| < £1. 

Therefore, for all t > T 1; we have 

<x(£) <x*(t)+ei. (36) 
Thus, by the second equation in (flj and by Bj4j we conclude that 
I'(t) <p(t,S(t),R(t),Q(t),I(t)) 



I(t) ~ 1(f) 



(r/(t) + a(t) + d(t) + on(t)) 



and finally 



^ limsup <p(t,S(t) 0,0,6) _ + + + 
^ limsup b$(t, S(t)) 

6^0 



l °ZTFF\= / 7TV dT < I limsup 6i(*,5(s))ds. 



Therefore, by HQ), flU and © we have 



/* lim sup bs(s,S(s)) ds f* lim sup fc^ (s,a;* (s)+£i ) 



7(t)=/(Ti)e 1 ^° ^/(Ti)e Jl «-» 

By p4|) we conclude that 

J(t) </(Ti)e 



Therefore lim J(i) = and this proves that the infectives go to extinction. 

t— >+oo 

Next, still assuming that R e (X) < for some A > 0, we let (S(t),R(t), Q(t),I(t)) 
be a solution of ([TJ) with non-negative initial conditions and (So(t), 0, Ro(t), Qo(t)) 
be a disease-free solution of ([l} with non-negative initial conditions. By [mj in 
Proposition [TJ we can assume that (S(t), I(t), R(t),Q(t)), (S (t), 0, R (t), Q (t)) G 
A ^ for t ^ T . Thus 

So(t),R (t),Q (t),S(t),R{t),Q(t) < K, 

for t Tq. By Theorem [2] we have lim I(t) = 0. Therefore given e > there 

exists T e ^ T such that 1(f) < e, S(t), R(t),Q(t) < K for t > T e . Thus 

A(t) - d(t)S ^ S"(f) > A(f) - d(t)S - eiV, 

for all t ^ T e . By comparison we get 

xi(t) > S(f) 5* x 2 (f) (37) 

where Xi(f) is the solution of (@| and x 2 (i) is the solution of (|20[) with £i = e and 
x t (T e ) = x 2 (T e ) = S(T e ). By|vj) in Proposition [5] we obtain, for all f ^ T u 

\x 2 (t)-xi(t)\^eND, (38) 

Since So(t) is a solution of (j3J), it is globally uniformly attractive, by[Iv]) in Propo- 
sition^ Thus, there is T 2 ^ T £ such that 

\ Xl (t) - S a (t)\ ^ e (39) 
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for all f > T 2 . By fl3Z|), J38J and ([39]) we conclude that 

^ a; 2 (t) ^ - eA^L> ^ 5 (t) - (1 + jV£>)e, (40) 
for all t ^ T 2 , and also, by ([37]) and ([39]). 

S(t)<5o(t) + s (41) 



for all t^T2- Since e > can be made arbitrarily small by taking t ^ T e , by (|40]l 
and ([4l]) we have 

lim |5(t)-Sf (t)| =0. (42) 

£— >+oo 

Let T 3 ^ T e be such that \I(t) - Io(t)\ < e for every t ^ T 3 and write = 
<3(i) — Qo(t). By the third equation in ([1]), we get 

-(d S + Ms + e s )u{t) < u'(t) < £(75 — (dj + (0.2)1 + Ei)u(t). 

We obtain 

u(t) > u (T 3 )e" (ds+(Q2)s+es)( *~ T3) 

and 

u(t) ^ u(T 3 ) e -(di+Mr+ei)(t-T 3 ) + ^ s (i e j _ 

or + (0-2)1 + £/ 

and thus 

lim |Q(f)-Q„(*)| = Hm u(t) < *° s . (43) 

t->+oo t->+oo dj + (02)1 + £l 

Since e > can be made arbitrarily small we get 

lim \Q(t)-Q (t)\ =0. 

t— >-+oo 

Let T4 ^ Ti be such that \Q(t) — Qo(t)\ < £ for every t ^ T4 and write w(t) = 
i?(t) — Ro(t). By the fourth equation in ([T]) and since 2o(i) = and for all t ^ 0, 
we get 

-d s w(i) s; w'(t) s: e(j s + e s ) - d lW (t). 

We obtain 

w(T A ) e-^Ct-TO <; wW <; w(T4) e - d x(t-T 4 ) + ^75 + £ S ) (1 - ^0 ) 
and thus 

hm |i?(<) - flo(t)| = lim «,(t) < £(75 . + £g) . (44) 

Since e > can be made arbitrarily small we get 

lim \R(t) - R (t)\ = 0. 
t— > +00 

By (|42f , ([43]) , (|44| and since all the infectives go to extinction, we conclude that 
any disease-free solution (So(t), 0, Ro(t), Qo(t)) with nonnegative initial conditions 
is globally attractive. We obtain [2] in Theorem [1] and this concludes our proof. 
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9. Simulation 

We carried out some experiments to illustrate our results. We considered some 
nonautonomous and nonperiodic models with mass-action and quarantine adjusted 
incidences. 

Firstly consider a nonautonomous model with mass-action and parameters A = 
0.001, d = 0.035, 7(t) = 0.4, a(t) = 0.01, ai(i) = a 2 (t) = s(t) = 0.2 and /3(i) = 
a(l - 0.7sin(0.3i))(2 - e -t ). For any A > 0, we obtain 

i?f 7M (A) =liminf 0.0886a(A - 4.667 sin(.3i + .15A) sin(.15A)) 

and 

i?f /M (A) = limsup 0.0886a(A - 4.667 sin(.3i + .15A) sin(.15A)) 

t— too 

(where the numbers R^ IM (A) and Rf IM {\) are given in section [7]). For instance, 
setting a = 9 and A = 21, we have R% IM {\) = 1.10599 > 1 and we can con- 
clude that the infectives are permanent and, setting a — 8 and A = 21, we have 
Rf IM r (A) = 0.98310 < 1 and we can conclude that the infectives go to extinction. 
Figures Q] and [2] illustrate the referred situations. 

0.0030 p 7.x 10-' r 




Figure 1. R S P IM = 1.10599 Figure 2. R s e 1M = 0.98310 



Now consider a nonautonomous model with quarantine adjusted incidence and 
the same parameters: A = 0.001, d = 0.035, j(t) = 0.4, a(t) = 0.01, a x {t) = a 2 (t) = 
e(t) = 0.2 and /3(t) = a(l - 0.7 sin(0.3t))(2 - e - '). For any A > 0, we obtain in this 
case 

R>? /QA (\) = liminf 3.10078a(A - 4.667 sin(.3< + .15A) sin(.15A)) 

and 

Re /QA W = limsup 3.10078a(A- 4.667 sin(.3t + .15A)sin(.15A)) 

t— too 

(where the numbers Rp^(X) and Rp^(X) are given in section[7]). Setting a — 0.25 
and A = 0.2, we have Rp /QA (X) = 1.07527 > 1 and we can conclude that the 
infectives are permanent and, setting a = .23 and A = 0.2, we have i?f /M (A) = 
0.989247 < 1 and we can conclude that the infectives go to extinction. Figures [3] 
and H] illustrate the referred situations. 
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